Draft version December 24, 2009 

Preprint typeset using L^TgX style emulateapj v. 08/22/09 



MORPHOLOGIES OF LOCAL LYMAN BREAK GALAXY ANALOGS IL A COMPARISON WITH GALAXIES 
AT Z :^ 2 - 4 IN ACS AND WFC3 IMAGES OF THE HUBBLE ULTRA DEEP FIELD 

R. A. Overzier\ T. M. Heckman^ D. Schiminovich^, A. Basu-Zych^, T. GoNgALVES^ D. C. Martin^ & R. M. Rich^ 

Draft version December 24, 2009 



On 
O 
O 
(N 

O 

Q 

o 

u 

6 

<N 

> 



o 



ABSTRACT 

Previous work has shown that Lyman break galaxies (LBGs) display a range in structures (from single 
and compact to more clumpy and extended) that is different from typical local star- forming galaxies. 
Recently, we have introduced a sample of rare, nearby {z < 0.3) starburst galaxies that appear to 
be good analogs of LBGs. These "Lyman Break Analogs" (LBAs) provide an excellent training set 
for understanding starbursts at different redshifts. We present an application of this by comparing 
the rest-frame UV and optical morphologies of 30 LBAs with those of galaxies at z ~ 2 — 4 in the 
Hubble Ultra Deep Field. We compare LBAs with star- forming sBzK galaxies at z ~ 2, and LBGs at 
z ~ 3 — 4 at the same intrinsic UV luminosity {Luv ^ 0-3i^*^3). The UV/optical colors and sizes of 
LBAs and LBGs are very similar, while the BzK galaxies are somewhat redder and larger. LBAs lie 
along a mass-metallicity relation that is offset from that of typical local galaxies, but similar to that 
seen at z ~ 2. There is significant overlap between the morphologies (G, C, A and M20) of the local 
and high redshift samples, although the high redshift samples are somewhat less concentrated and 
clumpier than the LBAs. Based on their highly asymmetric morphologies, we find that in the majority 
of LBAs the starbursts appear to be triggered by interactions/mergers. When the images of the LBAs 
are degraded to the same sensitivity and linear resolution as the images of LBGs and BzK galaxies, we 
find that these relatively faint asymmetric features are no longer detectable. This effect is particularly 
severe in the rest-frame ultraviolet. It has been suggested that high redshift galaxies experience intense 
bursts unlike anything seen in the local universe, possibly due to cold flows and instabilities. In part, 
this is based on the fact that the majority (~70%) of LBGs do not show morphological signatures of 
interactions or mergers. Our results suggest that this evidence is insufficient, since a large fraction 
of such signatures would likely have been missed in current observations of galaxies at z ~ 2 — 4. 
This leaves open the possibility that clumpy accretion and mergers remain important in driving the 
evolution of these starbursts, together with rapid gas accretion through other means. 

Subject headings: cosmology: observations - early universe - galaxies: high-redshift - galaxies: star- 
burst 



1. INTRODUCTION 

One of the key tasks in galaxy evolution is to un- 
derstand how the young, forming galaxies observed 
at high redshift relate to the well-defined Hubble se- 
quence observed at the present epoch. The study 
of sizes and morphologies of large samples of galaxies 
as a function of redshift would not have been possi- 
ble without the Hubble Space Telescope (HST). Stud- 
ies of the rest-frame UV sizes of Lyman Break Galax- 
ies (LBGs) at z ~ 3 — 7 indicate that they are 
mostly very compact objects with a single core (ri/2 ~ 
0.7 — 1.5 kpc), and that the size distribution devel- 
ops a tail of larger sized objects of up to several 
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kpc towards lower redshift s (e.g. [B ouwens et al.' '2 0041: 
iFerguson et al. 2004; Papo vich et a l. 2005; Oesch e t al.l 
l2009b) . Morphological studies performed by means 
of (a combination of) visual classifications, quantita- 
tive morphologi cal parameters and two-dimensional pro- 
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file fi tting (e.g. "Abraham et al.' '1996"; 'Lotz et al 
|2006l: [E lmegreen et al. 2005; Ravindranath et al 
Consel ice fc Arnoldii2009l : iPettv et al.n 2009) have shown 
that this increase in sizes at 2; ^ 2 — 4 is related to the 
accumulation of luminous star-forming clumps or cores 
within a variety of structures, including spheroid- and 
disk-like objects and irregular objects. These clumpy 
systems are expected to coalesce and form a spheroid 
while a surrounding disk may grow throu gh the con- 
tinued accretion of gas (jElmegreen et al.l l2008l . l2009l : 
iDekel et al.ll2009l ) . Individual clumps could in some cases 
be the nuclei of star- forming objects that are merging to- 
gether, or they could be giant starburst regions inside a 
larger gaseous system induced by merging or a plenti- 
ful smooth or "lumpy" gas supply. Detailed knowledge 
on the importance of such processes would, in principle, 
provide powerful constraints on models of galaxy forma- 
tion ( e.g. iBirnboim fc De kel"2003 |:rSomerville et al.l l2001l. 
[2OO8I : IGuo White! 12008. 2009r iDekel et al.ll2009[ ). but 
observationally they are hard to ascertain, especially at 
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high redshift. First, estimates of the (major) merger rate 
of galaxies by means of galaxy pair counts are difficult 
and critically depend on the mer ger time-scale, which 
may evolve with redshift (Kitzbi chler fc White! [2008). 
Second, our ability to identify galaxy mergers is a strong 
function of, e.g., the stage o f the merg er, the viewing 
angle, and gas fraction (e.g. iLotz et al] [2008). Third, 
at high redshift it is neither possible to directly mea- 
sure the (HI) gas fractions of galaxies nor to map the 
distribution of intergalactic hydrogen gas believed to be 
the supplying reservoir. However, the relatively wide 
range in kinematic properties of the emission line g:as 



observed in high redshift galaxies (e.g. [Law et al 



20091: iLehnert et all 120091 : iForster Schreiber et al 



200^, 



20091; 



Lowenthal et al.|[2009l ) may indicate that a variety of the 



basic mechanisms outlined above could be at play. 

One of the most basic tools that can be used to study 
the origin of these peculiar morphologies at high red- 
shift is to contrast them agai nst local or lower red- 
shift sa mples of galax i es (e.g. [Hibba rd fc Vacca' "1997'; 
Burgarell a et all 120061: iLotz et al.ll 2006; Scarlata et al. 
2007i : iCardamon e et all l2009i nElmegree n et al.i i2009l : 



Ostlin et a l. 2009i; iPettv et all [2009 : Rawat et al.[[2009[ ). 



However, it is important to keep in mind that the prop- 
erties of galaxies in the local universe are generally very 
different from those at high redshift, making it hard 
to disentangle actual physical differences from observa- 
tional biases. To better facilitate the straight compar- 
ison, the "Lyman break analogs" (LBA) project was 
designed in order to search for local starburst galax- 
ies tha t share typical characteristics of high redshift 
LBGs ([Heckman et al.[ [2005[ ). In brief, the UV imag- 
ing survey performed by the Galaxy Evolution Explorer 
(GALEX) was used in order to select the most luminous 
{Lpuv > 10-^^-^ Lq) and most compact {Ipuv > 10^ 
Lq kpc~^) star-forming galaxies at z < 0.3. Although 
such galaxies are very rare, they tend to be much more 
luminous in the UV than typical local starburst galaxies 
studi ed p reviously (Heckman et aLj 119981 : IMeu rer et al. 
[19971 . [1999 ) consistent with high SFRs and relatively lit- 
tle dust extinction. The median absolute UV magni- 
tude of the sample is -20.3, corresponding to 0.5L*^3, 
where is the cha racteristic luminosity of LBGs at 

z - 3 fSteidel e t all [19991 . i.e., Muoo^ab = -21.07). 
Analysis of their spectra from the Sloan Digital Sky Sur- 
vey (SDSS) and their spectral energy distributions from 
SDSS, GALEX, Spitzer and VLA and follow-up spec- 
troscopy subsequently showed that the LBAs are similar 
to LBGs in their basic global properties, including: stel- 
lar mass, metallicity, dust extinction, SFR, and emission 
line gas properties ([H eckman et al.[ [20 05: Hoo pes et all 
[20071 : [Basu-Zvch et ai.L.2007: Overzier et al. 200i r 

One of the advantages of LBAs being so bright in the 
rest-frame UV is that their morphologies can be easily 
compared with typical galaxies at high redshift without 
having to artifically brighten th em as was d one in pre- 
vious studies (Lotz efal] [20041 : [Pettv et all l20Q9l). In 
[Overzier et al.[ (|2008[ ) (Paper I) we analyzed the UV 
morphologies of a small sample of 8 LBAs observed 
with HST, finding that most of the UV emission in the 
LBAs originates in highly compact burst regions in small, 
clumpy galaxies that appear morphologically similar to 
LBGs. We argued that if LBGs at high redshift are also 



small merging galaxies similar to the LBAs, this would 
be very hard to detect given the much poorer physical 
resolution and sensitivity. In this paper, we present an 
analysis of the morphologies of our full data set consist- 
ing of rest-frame UV and optical HST images of 30 LBAs 
and compare with star- forming galaxies at z c:^ 2,3,4 
in the Hubble Ultra Deep Field (HUDF). The structure 
of this paper is as follows. In Section 2 we describe 
our low and high redshift samples, the observations and 
data reduction methods, and our techniques for produc- 
ing redshifted simulated images as well as for perform- 
ing parametrized galaxy morphologies. In Section 3 we 
compare the rest-frame UV and optical colors, sizes and 
morphologies of the LBAs, BzKs and LBGs. We discuss 
our results in Section 4, followed by a summary of the 
main results. We use the AB magnitude system through- 
out the paper, and assume a cosmology [Qm^^a^Hq] = 
[0.27,0.73,73.0] (with Hq in km s"^ Mpc"^) so that the 
angular scales at z ~ 0.2 and 3.0 are about 3 and 8 kpc 
arcsec"^, respectively. 

2. SAMPLES, DATA, AND SIMULATIONS 
2.1. Lyman Break Analogs at z < 0.3 

We use the data from our sample of 30 LBAs ob- 
served with HST (Programs 10920/11107). Rest-frame 
UV images for 24 objects were obtained with the Ad- 
vanced Camera for Surveys (ACS) Solar Blind Channel 
(SBC) in the filter F150LP (A^ ^ 1614A, twenty-four 
objects) and with the High Resolution Camera (HRC) 
in the filter F330W (A^ ^ 3334A, six objects). The ex- 
posure times per object were typically ~2500 s. Match- 
ing observations in the rest-frame optical were carried 
out either using the Wide Field and Planetary Camera 
2 for -3600 s in F606W (A^ ^ 6001A; twenty-four ob- 
jects) or using the ACS Wide Field Channel for -2200 
s in F850LP (A^ ^ 9170A; seven objects). The ACS 
and WFPC2 observations were divided into, respectively, 
three and six dithered exposure s and were combined us- 
ing multidrizzle (Koekemoe r et al.[[20Q2f). These data 
were presented in Overzier et al.[ (|2008[ . [20Q9IV The UV 
and optical images are shown in Fig. [H 

2.2. Star-forming Galaxies z :^ 2,3,4 

In order to be able to compare the morphologies of 
LBAs with the best quality morphological data possi- 
ble, our high redshift samples all come from the Hubble 
Ultra Deep Field (HUDF). The comparison data con- 
sist of three samples of star- forming galaxies at 2: — 2, 
z — 3 and 2: — 4. These three samples were chosen 
because starburst galaxies can be efficiently selected at 
each of these redshifts, using photometric color criteria 
described below. This has resulted in good statistical 
data on the sizes, morphologies and stellar populations 
of such galaxies. 

Our highest redshift samples are based on the pho- 
to metric r edshift catalog of galaxies in the HUDF from 
iCoe et al.[ (j2006). We selected two samples of LBGs hav- 
ing redshifts in the range 2.5 < z < 3.5 {z — 3) and 
3.5 < z < 4.5 (z — 4). We only selected objects having a 
high photometric redshift accuracy ("ODDS" >0. 99, i.e., 
the Bayesian redshift probability distribution, P{z), is 
peaked around the adopted best-fit photometric redshift 
without any significant secondary peaks or plateaus at 
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Fig. 1. — HST images of LBAs in the UV and optical at their true and similated redshifts. For each object, we show a set of four postage 
stamps. The top-left and top-right show the observed UV (6^^ x 6^^) and optical (10^^ x 10^^) images, respectively. The ID numbers and 
redshifts from Overzier et al. (2009) are indicated in the optical ('OPT') stamps. Below these images we show the same object simulated 
at 2 = 3 and observed at the HUDF depth through the ACS/WFC Veoe and W^FC3/IR ifieo filters probing the rest-frame UV (bottom- 
left) and optical (bottom-right), respectively. The simulated images measure 2'' x The observed UV images most often show one to 
several high surface brightness clumps, while in many of the optical images faint tidal features and companions suggest that the UV-bright 
starburst regions are triggered by mergers between small galaxies. In most cases the UV and optical light distributions are very similar, 
although in a small number of cases (e.g., objects 080844 and 210358) the UV emission is dominated by a single point-like source situated 
in a more extended and often disturbed disk seen in the optical (these objects were referred to as Dominant Central Objects (DCOs) in 
lOverzier^t al. (2009)). In most cases, the loss in sensitivity and resolution in the artificially redshifted images causes the UV clumps to 
blend into larger regions and hides the underlying structures and subtle features associated with the mergers seen in the optical. For color 
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other redshifts). It is important to note that we are 
not working with spectroscopic redshifts, meaning that 
we anticipate some scatter between the actual and the 
adopted redshifts^. Nonetheless, after we have applied 
an additional cut on UV luminosity (see below), it is 
expected that these samples are very similar to the typ- 
ical UV-selected LBG samples referred to in literature 
as tZ-dropouts {z ~ 3) and 5-dropouts {z ~ 4). At 
z ~ 2, analogous selections can be made based on the 
UGR colors of starbursts that are typically referred to 
as 'BM/BX' samples. Because deep U and G-band data 
were not available to us for the HUDF, we have instead 
used a sample of galaxies selected accordin g to B zK = 
{z-K)-{B-z)> -0.2 from lKong eHIl f2QQ8h . This 
criterion selects star- forming galaxies ( ^^sBzK^^ galaxies) 
with photometric redshifts in the range z ~ 1.4 — 2.5 
([Daddi et al."2004), and it has been shown that there 
is a high degree of overlap between the UV-selected 
"BM/BX" galaxies and the i^-selected sBzK galaxies, 
particularly at faint (Kvegg ^ 22) magnitudes (e.g. 
Reddv et al]l2QQ5L I2QQ6I : iGrazian et al]l2QQ7l : iKong et al.l 



2QQ8f ). However, we note that the sBzK selection is much 



more complete in terms of dust extinction, and thus in- 
cludes a higher fraction of relatively obscured galaxies 
missed by the UV selections. In any case, in order to 
make sure we are comparing objects of similar intrinsic 
UV luminosities to those of the LBAs, we have limited 
the sBzK and LBG samples to Muv,i700 < —19.5 in our 
analysis below. After further cleaning the samples from 
objects that are blended or near bright foreground ob- 
jects or that lie too close to the image edges, we are left 
with 30 {z ~ 2), 66 {z ~ 3) and 45 {z ~ 4) objects in 
the region of the HUDF covered in the optical imaging 
provided by the ACS. 

To measure morphologies in the res t- frame far-UV 
we us e the ACS Veoe-band image from iBeckwith et al.l 
(|2QQ6f ). These data cover an area of 11.2 arcmin^ and 
were taken in 122 dithered exposures with a combined 
exposure time of 135320 s. The frames were drizzled 
at an output pixel scale of 0'.'03 pixel" ^ with a PSF of 
full-width-at-half-maximum (FWHM). In order to 
also measure morphologies in the rest-frame optical at 
high redshift we require high quality near-infrared data. 
We use data from the WFC3/IR channel observations 
performed as part of the early release science observa- 
tions program 11563 (Bouwens et al.l l2009al : lOe sch et al." 
^009a). We only use images taken through the filter 
^160 (Ac ^ 1-54 /im) covering an area of 4.7 arcmin^ 
centered on 3^32^38.5" and -27^4^0.0''. In order to 
create the i^ieo combined mosaic image, we started from 
the pipeline calibrated (fit) images released as part of 
the Servicing Mission 4 Early Public Observation Data 
Products as input to MultiDrizzle. The i^ieo data were 
taken in 14 visits of 2 orbits each during August 26, 2009 
to September 6, 2009. After removal of several visits 
and exposures of reduced quality, we were left with 47 



^ We further note that the lCoe"et al. (2006) photometric redshift 
analysis did not incorporate any [/-band data, which is important 
for a careful sampling of the Lyman break at z ~ 3. However, 
it has been shown that the main advantage of the inclusion of 
tZ-band data is that it resolves some ambiguities in the P{z) of 
genuine 2; ~ 3 objects, while there are relatively few low redshift 
interlopers that are wrongfully placed at z ^ 3 when ?7- band data 
is not available (jNonino et al.|[2009l : [Rafelski et al.ll2009V 



TABLE 1 

Overview of Rest-frame Wavelengths of the 
Observations/Simulations. 



LBAs (N") 


BzKs 


LBGs 


LBGs 


z < 0.3 


z ~ 2 


zr^3 


Zr^4 



F150LP (24) 
-I350A 



Rest-UV 

ACS/WFC F606W 
'2000A ~1500A ~1200A 



F606W (23) 
~5000A 



F330W (6) 
^2800 A 

Rest-optical 



F850LP (7) 
-75OOA 



WFC3/IR F160W 
-5IOOA ~3900A - 



^ Numbers between parentheses indicate the number of LBAs 
observed in each filter. 

dithered fit images with a total exposure time of ?^66000 
s. Two passes through multidrizzle were made. Dur- 
ing the first stage, we created single distortion-corrected, 
registered images and combined them into a median im- 
age. Because the median image contained significant lev- 
els of background structure not removed by the pipeline, 
SExtr actor was used to create an object-free background 
image. This background image was blotted back to the 
geometry of each of the original input frames and sub- 
tracted. In the second stage, the frames were drizzled to- 
gether to produce a final cleaned background-free mosaic 
with an output pixel scale of ff!06 pixel" ^ (pixfrac=0.7) 
and a PSF of - ff!2 (FWHM). Because of the smaller 
coverage in the near-infrared provided by the Wide Field 
Camera 3 (WFC3), the sample sizes are 22 at z ~ 2 and 
35 at z ~ 3. 

In Section [3] we will also make use of some of the pho- 
tometric data in the ACS filters ^435, Vbob^ ^77^ a - nd z^^ o 
taken from the HUDF catalog from iCoe et all (|20Q6( ). 
and groundbased H- and i^Tg-band data from the pub- 
licly available Great Observatories Origins Deep Survey 
Multiwavelength Southern I nfra red Catalog (GOODS- 
MUSIC) from lSantini et all (|2QQ9f ). 

2.3. Redshift Simulations 

In order to be able to compare the sizes and morpholo- 
gies of our low and high redshift samples, we performed 
redshift simulations of the \ow-z LB A sample. We apply 
corrections for cosmological surface brightness dimming 
and for changes in physical resolutions. The first step of 
the procedure is to rebin the low-z images by a factor 
b = (^1/^2) (52/51), where 6i is the angle on the sky of an 
object of fixed size d at z = z^, and Si is the instrumen- 
tal pixel scale (in arcsec pixel"^). The rebinning factor 
can be expressed in terms of the ratio of angular diame- 
ter distances, b = {D A2 / D Ai){s2 / si) ^ where = d/Oi 
and Z2 > zi. The second step is to reduce the surface 
brightness of each (rebinned) pixel according to the rel- 
ative amount of cosmological dimming of a galaxy at 
Z2 with respect to that at zi. We calculate the scaling 
by making use of the fact that the absolute rest-frame 
magnitude (or luminosity) of the object before and after 



redshifting will be conserved (M^ 



A2/(l + 22) 



M 



Ai/(l+2i) 



with matched filters so that A2 ~ Ai(l + Z2)/(l + 2^1)). 

The effectiveness of our simulations is limited mainly 
by two factors. One, for each LB A we have only a single 
band in each of the UV and optical available. This means 
that there will be uncertainties in the (morphological) K- 
corrections. While our rest-frame far-UV images are a 
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Fig. 2. — Panels show the distributions in absolute UV magnitude (left), UV continuum slope puy (middle) and rest-frame UV-optical 
(FUV — R) color (right). In each panel, black histograms correspond to the LBAs, while blue, green and red histograms show the 
distributions for the z ~ 2, z ~ 3, and z ~ 4 comparison samples, respectively. The vertical dotted dashed line in the left panel indicates 
the magnitude limit applied to the comparison samples. On the top axis of the middle panel we give the inferred E(B — V) extinctions 
corresponding to the measured ^5 of a young star-forming galaxy (Calzetti 2001). On the top axis of the panel on the right we give the 
inferred age corresponding to the color measured for a young star-forming galaxy of solar metalliticy and a reddening of E(B — V) = 0.1 
observed at z = 0.2. When calculating the UV and UV-optical colors, we did not perform filter interpolations but instead for each sample 
we used those filters that were closest in rest-frame wavelength to the desired colors. The four samples are largely comparable in all three 
properties, although the median UV and UV-optical colors of the BzK sample is somewhat redder compared to those of the LBAs and 
LBGs. The histograms in the middle and right panels were all normalized to the area of the LBA histograms. 



near perfect match to the observed optical ACS images 
of the LBGs, the six [/-band images probe rest-frame 
wavelengths that are slightly redder even than the ACS 
Z850-band data on LBGs (WFC3/IR F105W and F125W 
would provide a better rest-frame match). Second, our 
rest-frame optical images probe rest-frame wavelengths 
of ~ 5000 — 8000A, while the reddest band available 
for LBGs in the HUDF (WFC3/IR F160W) only probes 
the rest-frame optical at A ~ 4000A even for LBGs at 
z ^ 3. We will remedy these shortcomings as follows. 
In the minority of cases where our available UV data 
is somewhat bluer or redder than the mean rest-frame 
wavelenghts of the high redshift data, we will assume 
that the morphological i^-corrections from rest-frame 
FUV to NUV/[/ are negligible. This is a reasonable 
assumption given that the morphological i^-corrections 
for both LBAs and LBGs are found to be quite modest 
even from the rest-frame UV to the rest-frame optical 
(see, e.^.. iLotz et al.l [2QQ1 I2QQ6I : Papovich et al. 2005; 
iConselice fc Arnoldl 12009. and Section [3] of this paper). 
We will further assume that our relatively red rest-frame 
optical images of LBAs are always representative for the 
morphologies at A > 4000A. We then simulate rest-frame 
UV images in F606W (ACS/WFC) and rest-frame opti- 
cal images in F160W (WFC3/IR), both at the HUDF 
depth. However, in all cases we apply appropriate color 
terms based on the full UV-optical spectral energy dis- 
tribution that is available to us from GALEX and SDSS 
in order to minimize at least the spectral i^-corrections. 
This will ensure that the objects have the correct surface 
brightness when artificially redshifted. Because of the 
relatively blue rest-frame wavelength of F160W we will 
only simulate rest-frame optical images for comparison^ 

^ In a future paper we plan to perform similar simulations for 
JWST or AO-assisted, ground-based observations in the observed 
K-band in order to compare morphologies of LBAs and LBGs at 
longer rest-frame optical wavelengths. 



with LBGs at z ~ 3. Our simulated images take into ac- 
count the effects of sky background, Poisson noise, dark 
current, readnoise, sub-exposures, and PSF convolution 
based on stars in the actual HUDF Vgog and i^ieo im- 
ages. The different instrumental configurations and cor- 
responding rest-frame wavelengths used for the observed 
and simulated data are summarized in Table [H Exam- 
ples of the ACS/WFC Vqog and WFC3/IR i^ieo images 
simulated based on LBAs redshifted to z = 3 are shown 
in Fig. [TJ 



2.4. Size and Morphology Measurements 

We will compare the galaxy radii r^o and rgo contain- 
ing, respectively, 50 and 90% of the light measured in 
the UV and optical using SExtractor. The "total" light 
radius was set to 4 times the Kron radius. We will 
also compare a number of morphologial quantities used 
widely in the literature: the Gini coefficient {G; a mea- 
sure of the equality with which the flux is distributed 
across a galaxy), M20 (the log of the ratio of the second 
order moment of the pixels containing the 20% bright- 
est flux to the total second order moment), concentra- 
tion (C; five times the log of the ratio of the circular 
radii containing 80 and 20% of the flux), and asymme- 
try {A; a measure of the mirror symmetry of an object). 
A fifth and commonly used parameter that measures the 
dumpiness {S) of galaxies will not be used here, as it was 
found to be of limited use at high redshift where galax- 
ies te nd to be faint and compact (Conselice fc Arnoldl 
I2OO9I). We also do not perform any two-dimensional pro- 
file fitting, but note that inferences made based on such 
studie s applied to LBGs can be found elsewhere (e.g. 
Ravin dranath et al.l[2006l : iRawat et aLll2009l : iPettv et al.l 
2009). To calculate the four parameters we closely fo llow 
the definitions and procedures described in .Lotz et al.l 
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(j2QQ4 I2QQ6D and Paper I and in the footnote^. Some 
LBAs have faint companions that were included in the 
simulated images. For the high redshift samples, it is 
often hard to determine whether neighbouring objects 
are physically associated or not given the general clumpy 
nature of BzK and LBG galaxies as well as overcrowd- 
ing in the deep HUDF images. Our SExtractor settings 
were chosen such that high surface brightness regions 
connected by diffuse emission were largely considered as 
a single system, while other neighbouring objects were 
masked out fairly aggressively. This should be kept in 
mind as the inclusion and rejection of close neighbours 
can have significant effects on some of the morphological 
parameters. However, we believe that our main conclu- 
sions will not be affected by this. 

3. RESULTS 

3.1. Rest- frame UV- optical Colors 

In order to re-emphasize the high degree of similarity 
between the LBAs and typical star-forming galaxies at 
high redshift, we compare the main UV and optical pho- 
tometric properties in Fig. [2l In the panel on the left, we 
first show the distribution in absolute UV magnitudes at 
1700 A for the LBAs (black histo gram) and the three high 
redshift samples {BzK in blue, [/-dropouts in green and 
5-dropouts in red). As stated in Section [2^ we match 
the four samples in UV luminosity by placing a cut at 
Muv = —19.5 mag. In the middle panel of Fig. [2]we plot 
the distributions of the UV continuum slope Puv^ where 
fx (X and p is typically measured for A in the range 
1650A to 2300A co rresponding to approximat ely the far- 
and near- UV (e.g. iMeurer et all IT995I . [19991 ). For the 
LBAs, we estimate /3 from the (FUV-NUV) color, while 
we use (^435-^606) at 2; ^ 2, (Veoe-^rrs) at z ~ 3 and 
(^775~^85o) at 2: ~ 4. The top axis of the middle panel 
shows how a given /3 corresponds to the approximate 
E{B — V) extinction for a young star- forming galaxy 
(~100 Myr old and forming stars at a constant rate) as- 
suming the dust recipes from Calzetti (2001). The distri- 
butions are all very blue, consistent with no or only small 
amounts of reddening due to dust. The median redden- 
ing in the BzK sample is somewhat higher than that 
i n the LB A /LBG samples consistent wi th other works 
(jKong et all 120081 : iBouwens etHI l2009b[ ). In the right 
panel of Fig. [2] we show the distributions in rest-frame 
UV-optical colors for the four samples. The colors of the 
high redshift samples were calculated using the filters 
that were closest to the rest-frame central wavelengths 
probed by the {FUV - R) color of the LBAs at z - 0.2. 
On the top axis of the panel on the right we give the 

^ In brief, we use SExtractor to make an object segmentation 
map and mask out neighboring objects. The image is background 
subtracted, and we calculate an initial Petrosian radius (rp with 
?7 = 0.2) using the object center and (elliptical) shape information 
from SExtractor. We then smooth the image by a = rp/5 and 
create a new segmentation map by selecting those pixels that have 
a surface brightness higher than the mean surface brightness at the 
Petrosian radius. We recalculate the object center by minimizing 
the second order moment of the flux, and then recalculate the 
Petrosian radius in the original image using this center. The total 
flux is deflned as the flux within a radius of 1.5 x rp. C is calculated 
in circular apertures containing 20 and 80% of the light. A is 
calculated within a circular region of radius 4xr5o, and we subtract 
the asymmetry of the background using a similar sized region free 
of objects. 
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Fig. 3. — The stellar mass-metallicity (M*-Z) relation. Con- 
tours indicate the distribution from a large sample of star-forming 
and composite objects at z < 0.3 extracted from the SDSS/DR7 
spectroscopic sample. Large black symbols indicate the values mea- 
sured for individual objects in our LB A sample from Overzier et al.l 
(2009). Blue symbols indicate the average values found for star- 
forming galaxies at z ~ 2 in several mass bins from Erb e t al.l 
(2006). The grey points indicate the M^-Z distribution of a sta- 
tistical sample of ~500 LBAs extracted from the GALEX/GR5 
and SDSS/DR7 crossmatched samples (Overzier et al. 2010, in 
prep.). Objects identifled as star-forming/ AGN composites have 
been indicated using a darker shade of grey. All metallicities have 
been estimated using the "N2" method based on the [NII]/Hck; ratio 
from Pettini Sz Pagel (2004). The horizontal dotted line indicates 
solar metallicity, and we note that the metallicities derived using 
the N2 method saturate near this line. 

inferred age corresponding to the color measured for a 
young star-forming galaxy of solar metalliticy and a red- 
dening of E{B -V)= 0.1 observed at z = 0.2. The 
distributions are similar to those expected for galaxies 
having ages that peak around a few hundred Myr to a 
Gyr. 

3.2. Mass- Metallicity Relation 

The (stellar) mass-metallicity (M* — Z) relation is 
a crucial diagnostic for evaluating the gas-phase metal 
abundance of galaxies as a fun c tion of their ba ryonic 
mass (e.g. IZaritskv et al.l 119941 : iTremonti etld] l'2004l : 
ISavadio et al.ll2005f ). In Fig. [3] we reproduce the - Z 
relation of LBAs first shown bv lHoopes et al.l (|2007l ). and 
updated according to our most recent samples and mea- 
surements. Large black symbols indicate the values mea- 
sured for the 30 LBAs from our sample and using the 
SDSS-based stellar masses and emi ssion line ratios based 
presented in lOverzier et al.l (|2009l ) . Contours show the 
density distribution from a large sample of star-forming 
and composite objects at z < 0.3 extracted from the 
SDSS/DR7 spectroscopic sample^^. Although we do not 

Available on the website of the Max-Planck Institute for As- 
trophysics: http : //vrvrvr . mpa-garching . mpg . de/SDSS/DR7/ 
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have the gas-phase metallicities of the three main com- 
parison samples at high redshift that are used in this 
paper, we can at least compare with the mass-binned 
averages found for star-forming galaxies at z ~ 2 from 
lErb et all ([2006) (shown in blue). This z - 2 sample 
is very similar to our z ~ 2 BzK sample in most as- 
pects. All metallicities in Fig. [3] have been estimated 
using the ^^N2" m ethod based on the [Nil] /Ha ratio from 
iPettini fc Pagell (2004). The horizontal dotted line indi- 
cates solar metallicity, and we note that the metallicities 
derived using the N2 method saturate near this line. 

Fig. [3] illustrates the very similar distributions both in 
stellar mass and in metallicity for LBAs and the z ^ 2 
sample, and in such a way that their — Z relation is 
increasingly offset from the local relation tow ards lower 
stellar masses. As shown bv lErb et al.l (j2006l ). the offset 
with respect to the local relation is similar at all masses, 
and could be explained by the (mass-independent) loss 
of metals from supernova winds. Except for the highest 
masses, the LBAs show very similar offsets with respect 
to local galaxies of the same mass. This could either 
indicate that LBAs are still in the process of convert- 
ing a relatively large gas mass into stars compared to 
other galaxies of the same mass, or that they have had a 
recent accretion event of metal-poor gas associated with 
the onset of the starburst, possibly coupled with a (mass- 
dependent) outflow o f metals due to winds ( iHoopes et al.l 
120071 : lOverzier et all [2009). 

3.3. Sizes 

In Fig. m we show the half-light radius distributions 
for the four samples. The three panels on the left 
marked 'UV show the rest-frame UV sizes for each of 
the high redshift samples, while the two panels on the 
right marked 'OPT' show the rest-frame optical sizes for 
the z ^ 2 and z ^ 3 samples. In each panel, the sizes 
measured for LBAs simulated at corresponding redshifts 
are shown in black, while their "true" size distribution 
measured from the full-resolution, low redshift HST im- 
ages is indicated by the grey lines. On the top axes we 
have indicated the corresponding angular scales. The 
size distributions of the four samples are comparable, 
with a notable exception being the, on average, larger 
sizes measured for the BzK galaxies in the rest-frame 
optical. Our results are consistent with earlier stud- 
ies of the size distributions at high redshift finding that 
star-forming galaxies at 2: ~ 2 — 4 are compact galax- 
ies having half-light radii in the range ~ 0.5 to a few 
kpc in the UV and ~ 0.5 to 6 kpc in the optical, consis- 
tent with more detai l ed analyses from the l iterature (e.g. 
iBouwens et al.""2004': 'Ferguson et al]l2004l: iTruiillo et al.l 
[2006; Dahlen et al. 2007; Oesc h et al.ll2009bD ~ 

The main purpose for comparing the main physical 
quantities summarized in Figs. [2 [3] and [4| that are 
not sensitive to redshift-dependent observational effects 
(i.e., luminosity, color, mass, metallicity and size) was to 
demonstrate the basic similarities between our samples. 
The results support our basic premise that the sample 
of LBAs is useful for investigating various other proper- 
ties of high redshift galaxies, such as their morphologies 
that are likely very sensitive to redshift effects and that 
can not be easily obtained from the available data. In the 
following subsection we will compare the morphologies of 
LBAs, BzKs and LBGs in order to investigate whether 



they are similar as well or perhaps notably different. 

3.4. Morphologies 

Our results on the morphologies (G, C, A and M20) 
measured in the rest-frame UV and optical images fol- 
lowing our methods outlined in Sect. 12.41 are presented 
in Figs. [5H3 and Tables [2] -[H Before we discuss our 
results, we will briefly describe what t he d ifferent pa- 
rameters ar e most sen sitive to (e.g., se e [Abraham et al.[ 
1996, 2003 :jLQtz_et al.[ l2004. 2006; Cons elice et al.[[2000[ : 
Consel ice fc Arnold [2Q09> : iLisker. |2008, and references 
therein for derivations and detailed discussions). C is 
a powerful diagnostic for separating relatively shallow 
(low C) versus highly concentrated (high C) light pro- 
files, such as late-type spirals and bulge-dominated galax- 
ies. Two issues can affect C. First, galaxies may contain 
unresolved components containing more than 20% of the 
light, thereby causing R20 to be overestimated and C to 
be underestimated. Second, a small offset between the 
estimated center of a galaxy and a highly concentrated 
nucleus will also cause C to be underestimated compared 
to its true value. Because both effects can only be coun- 
tered effectively for the redshifted LBA sample (for which 
sub-resolution and sub-sensitivity spatial information is 
available) but not for the real high redshift samples, we 
have not attempted to correct for these effects. The Gini 
coefficient is known to correlate with C, but differs in 
the sense that it is sensitive to any co ncentrated light 
independent o f its spatial distribution ([Abraham et al.[ 
[20031 : [Lotz et al.l i2004). Similar to C and G, the M20 pa- 
rameter is also sensitive to concentrated light. However, 
because every pixel is multiplied by the square of its dis- 
tance relative to the centre (a free parameter), it is more 
sensitive to spatial fiuctuations such as bright clumps or 
double nuclei from merging systems (high M20). The 
deviation from rotational symmetry, A, can distinguish 
between symmetric profiles such as those of early-type 
galaxies (low A) and those that are less symmetric such 
as spirals, irregulars and mergers (high A). 

Below, we will first investigate the effect of redshift 
on the measurement of morphologies in the LBA sample 
itself (Sect. I3.4.ip . Then, we will compare the redshifted 
LBA samples with the high redshift comparison samples 
(Sect. [3X2D . 

3.4.1. Investigating redshift effects on LBA morphologies 

In Fig. [5] we compare the values calculated for 
LBAs at their intrinsic redshifts {z ~ 0.2) with those 
calculated after redshift ing them to z = 2,3,4. Black 
points indicate the individual measurements at low 
redshift, while green points indicate the measurements 
at z = 3. Large crosses indicate the median values at 
z ~ 0.2 (black), z = 2 (blue), z = 3 (green) and z = A 
(red), and the length of the bars indicate the 15-85 per- 
centile ranges. We summarize the main results as follows: 

1. In both the UV (left) and optical (right), the median 
values of G, M20 and A decrease from low to high red- 
shift, with A being the most sensitive to redshift: both 
its median value and its scatter are significantly smaller 
compared to low redshift. The C does not change in the 
UV. The drop in the median values can be explained due 
to a combination of the loss in resolution and sensitivity. 



8 



R.A. Overzier et al. 



Rgo (arcsec) R50 (arcsec) Rgg (arcsec) 

0.0 0.1 0.2 0.3 0.4 0.50.0 0.1 0.2 0.3 0.4 0.5 0.0 0.10.2 0.3 0.4 0.5 0.6 





12 




10 




8 








6 


1 


4 




2 









I 



uv 

z~3 



uv 

z~4 



012340123401234 
R50 (kpc) R50 (kpc) Rgo (kpc) 



B 
I 



12 
10 
8 
6 
4 

2 




RgQ (arcsec) R5Q (arcsec) 

0.0 0.2 0.4 0.6 0.8 1.00.0 0.2 0.4 0.6 0.8 1.0 1.2 
21 ■ ■ ' ■ ^ ■ ■ ' ■ ■ ^ 







OPT 


1 OPT : 






z ~ 2 :,r 





















2 4 6 
R50 (kpc) 



8 



2 4 6 
R50 (kpc) 



Fig. 4. — Panels on the left show the physical half-light radii measured in the rest-frame UV images for BzK galaxies (blue histogram), 
LBGs at z ~ 3 (green histogram) and LBGs at 2; ~ 4 (red histogram). The black histogram in each panel indicates the size distribution 
of the LB As as measured from the images simulated at each of the redshifts of the comparison samples. The grey line shows the "true" 
size distribution of the LBAs measured from the full-resolution, low redshift HST images. The two panels on the right show the size 
measurements based on the rest-frame optical data in a similar fashion. The rest-frame optical sizes are typically twice as large as the 
rest-frame UV sizes. In the rest-frame optical, BzK galaxies are somewhat larger than both LBAs and LBGs. The top axes of each panel 
indicate the corresponding angular sizes. The sizes were measured using SExtractor and have not been corrected for the size of the PSF 
{FWHM ^ O'/l in the UV and ^ 0^2 in the optical). 
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Fig. 5. — Distributions of the morphological parameters G, C, A M20 for the LBAs at their true and simulated redshifts in the rest-frame 
UV (left) and optical (right). Black crosses indicate the 15-85 percentile ranges measured off the true (low-redshift) images. For comparison 
we indicate the ranges found for LBAs redshifted to z = 2 (blue crosses), z = 3 (green crosses), and z = 4 (red crosses). The z ~ 4 sample 
is not shown on the right as we lack a sufficiently red band to probe the rest-frame optical light. Small symbols indicate the individual 
measurements for the LBAs at their intrinsic (low) redshift (black points) and LBAs simulated at 2; = 3 (green points). 



The loss of faint structures causes the 80% hght radius to 
be underestimated (lowering C) and the flux to be more 
evenly distributed over the detected pixels (lowering G). 
Resolution causes bright, clearly separated clumps seen 
in the unredshifted images to blend (see Fig. [1]), thereby 
lowering M20 and A. 

2. If we focus on the M20 — A plane (lower- left panels) 
- which as we will see below, in principle, serves as the 
best diagnostic for identifying highly disturbed galaxies - 
we flnd that the relatively high asymmetries measured for 
the LBAs at low redshift are less affected by the redshift- 
ing in the optical image compared to the UV. Previous 
studies of morphologies have suggested criterions of, e.g, 
M20 ^ —1.1 and/or A > 0.35 as indicators for double 



nucleated (e.^-, mer^n^) or generally disturbed galaxies 
(|Lotz et all I2QQ6I : IConselice fc Arnold! [2009) . Based on 
the M20 criterion, very few LBAs fall in this category 
even at z ~ 0.2 (respectively 3 and 5 in the UV and op- 
tical), while none would be selected at z = 3. Based on 
asymmetry, most LBAs are classifled as disturbed both 
in the UV and optical at low redshift and in the optical 
at high redshift, while only a very small fraction would 
be classifled as such in the UV at high redshift. 

3.4.2. Comparing LBAs, BzKs and LBGs 

In Fig. [6] we present our measurements for the high 
redshift samples, while in Fig. [71 we compare the 
parameter distributions measured for BzKs and LBGs 
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Fig. 6. — Distributions of the morphological parameters measured in the rest-frame UV (left) and optical (right) for the high-redshift 
HUDF samples. The z ~ 2 BzK sample is shown in blue, the 2; ~ 3 LEG sample in green, and the 2; ~ 4 LEG sample in red. The 
z ~ 4 sample is not shown on the right as the WFC3 i^ieo filter is too blue for probing the rest-frame optical at A > 4000A. Dashed 
lines indicate standard criteria used to separate mergers or disturbed galaxies (e.g., A > 0.35 V M20 > —1.1 V G > — O.I4M20 + 0-33 V 
G > — 0.4A + 0.68, where V is the logical OR operator) from galaxies having more symmetric profiles taken from Lotz et al.l (|2004l . [20081 ) 
and iConselice Arnold (2009). 



with the LBAs simulated at the same redshift (barred 
squared and circles indicate the median and 15-85 
percentile ranges for comparison samples and redshifted 
LBAs, respectively). We observe the following trends: 

1. At each redshift, there are systematic differences 
between the median values of LBAs on one hand, and 
BzKs /LBGs on the other. However, this difference is 
often comparable or smaller than the (large) scatter on 
the mean. The differences as well as the scatters are 
smaller in the optical compared to the UV. 

2. BzKs /LBGs (large squares in Fig. [7j) are some- 
what less concentrated, have lower G and higher M20 
compared to the LBAs simulated at the same redshift. 

3. In the rest-frame UV, the BzK galaxies (blue 
squares in Fig. [7j) have the most notable offsets com- 
pared to the LBAs (simulated at z = 2, blue circles), 
in the sense that they are much less concentrated, have 
higher asymmetries and larger M20. In this respect, they 
also differ from :^ 3 — 4 LBGs (although part of this is 
likely a redshift effect). In the optical, the offsets are con- 
siderably smaller but show the same general trend that 
the difference between BzK galaxies and LBAs is larger 
than that between the z 3 LBGs and LBAs. 

4. In the UV, a significant fraction of sources in each 
of the high redshift samples have M20 > —1.1, compared 
to none of the sources in the redshifted LBA samples. 
In the optical, the M20 is also higher for the true high 
redshift samples compared to the LBAs, although it 
must be noted that the 15-85 percentiles lie entirely at 
M20 < —1.1 for all samples. 

In Tables [2]-|4] we have summarized the main re- 
sults presented graphically in the three figures discussed 



above. Table [2] lists the median parameters found for 
each sample, and we estimate a scatter on the mean from 
the sample variances. Tables [3] and [H gives a measure 
of the statistical significance of the difference in median 
values between each pair of samples. We define the ratio 

Rp = (P2 — ^1 ) / ^ cTp^ + o-p^ , where Pi and crp. are the 

median and standard deviation in the morphological pa- 
rameter P measured for sample i. A value of \Rp\ ~ 1 
then indicates that the difference is comparable to the 
scatter, while \Rp\ ~ indicates that the difference is 
small. For each parameter and combination of samples, 
we also give the level of significance that the null hy- 
pothesis that the two samples are drawn from the same 
parent distribution is true, as calculated from a two-sided 
K-S test. The statistical results are consistent with our 
qualitative conclusions summarized above. 

4. DISCUSSION 

4.1. Morphologies of LBAs at different redshifts 

In previous works we have shown that LBAs are 
clumpy starburst galaxies with peculiar morphologies 
that are most consistent with merging (Overzier et al.l 
I2QQ8I . 120091). Although we currently do not know the 
relative contribution from major and minor mergers, 
the diversity in LBA morphologies at least suggests a 
range in merger conditions and stages albeit that they 
all have in common that their starbursts are all very 
young (ag es of a few ten s of Myr), compact and UV- 
luminous ( Overzier et al.| [2Q09). We have shown that, 
both in the UV and optical, asymmetry appears to be a 
better indicator for this merger activity than M20. The 
medi an A is larger than t h e mer ger criterion of A > 0.35 
from Conselice fc Arnold! (|2009l ). while the median M20 
is much lower than the merger criterion of M20 ^ — 1.1 



10 



R.A. Overzier et al. 



0.70 
0.6S 
0.60 
0.5S 
0.50 
0.4S 
0.40 



1r 





uv c 

• LB As @z = 2.0 

• LB As @ z = 3.0 

• LB As @z = 4.0 

□BzKs @ z ~ 2.0 
□LBGs @ z ~ 3.0 
□LBGs @ z ~ 4.0 G 



I rfH 



I EJ^ 



- -I |Jr I 



-1.0 -1.5 -2.0 -2.5 



4 0.40 0.45 0.50 0.55 0.60 0.65 0.70 



0.50 
0.45 



it ■ ■ 




UJlC 
rpiipi 



OPT 

• LB As @ z = 2.0 

• LB As @ z = 3.0 



□BzKs @ z ~ 2.0 
□LBGs @ z ~ 3.0 



iL ill 



0.40 0.45 0.50 0.55 0.60 0.65 



Fig. 7. — Comparison of the rest-frame UV (left) and optical (right) parameters for (redshifted) LBAs and the high-redshift comparison 
samples. Barred squares indicate the median and 15-85 percentile ranges measured for the 2; ~ 2 BzK sample (blue), the z ~ 3 LBG 
sample (green) and the 2; ~ 4 sample (red), while barred circles of corresponding colors indicate the values measured for LBAs simulated at 
the corresponding redshifts. In the UV, the distributions show significant overlap in morphological parameter space, although the median 
G, C and M20 of the high redshift samples indicate that they are somewhat less concentrated and clumpier compared to the (artificially 
redshifted) LBAs. The strongest offsets are observed for the BzK sample. In the rest-frame optical, there is better agreement between the 
median values measured for the BzK, LBG and LBA samples compared to those measured in the rest-frame UV. 



(e.g. Lotz et aL"2QM lConselice fc A rnold"2009\ We red- 
shifted the LBAs to z = 2—4 and found a significant drop 
in the UV asymmetries such that the median A fell well 
below the merger criterion. In the optical, the drop was 
smaller and less significant. M20 became even smaller 
as the different star-forming clumps blended into single 
light concentrations at high redshift. We thus conclude 
that if galaxies similar to LBAs were observed in the UV 
at high redshift, all but a few would be mistaken as being 
relatively smooth and symmetric rather than identified 
as mergers. In the optical, we estimate that ~50% would 
b e iden tified as symmetric. 

iLotz et al.l ()2QQ8l ) investigated in detail the compli- 
cated relation between the ability to morphologically 
classify merging galaxies on one hand, and the physi- 
cal details of the merger on the other. As expected, 
the merger observability time-scale (i.e., the time dur- 
ing which a merger is identified as such compared to the 
total merging time) critically depends on, e.g., orbital 
parameters, gas fraction, dust, SN feedback, and view- 
ing angle. These results explain, for example, why the 
G, C, ^4, M20 system is much more effective in identify- 
ing the roughly two-thirds of local ULIRGs having more 
than one nuclei that are seen either just prior to final co- 
alescence or in pro jection a t first passage, than single nu- 
cleus ULIRGs (jLotz et al.ir2QQ4 12QQ8). In a si milar fash- 
i on, th is "merger observability" as defined bv lLotz et al.l 
()2QQ8l ) can also explain why a significant fraction of the 
LBAs do not appear to be particularly disturbed for some 
or all of the morphological parameters. However, in addi- 
tion to these effects, our simple simulations clearly show 
that for LBA-like galaxies observed at high redshift the 
decreased physical resolution and sensitivity are of equal 
or even greater importance than the physical details of 



the merger. 

4.2. Morphologies of LBAs, BzKs and LBGs 

The median M20 of LBGs at z ~ 3 — 4 is larger than 
that of LBAs, while the A is similar. The fractions of 
LBGs classified as disturbed based on either one of these 
parame ters in the UV is ^30%, consistent with previou s 
works ^Lotz et al.ll2QQ4 12QQ6I : IConselice fc Arnoldl[2QQ9[ ) . 
In the optical, a similar result is obtained based on A, 
while a much smaller fraction is derived based on M20. 
It is interesting to note that the differences observed be- 
tween BzK galaxies and LBAs are much more significant 
than those between LBGs and LBAs: in the UV, BzKs 
have the highest A and M20 and the lowest C and G of 
all samples considered here. In the optical, BzKs had 
the largest median half-light radii. 

It must be noted, however, that our way of selecting 
LBAs based on a UV surface brightness criterion was 
tuned to select objects having similar UV luminosities 
and sizes as typical LBGs at z ~ 3 (iHeckman et al.l 
l2QQ5f ). In principle, the LBA selection could be expanded 
or modified to search for such objects more similar to 
BzK galaxies, but one must be more careful as such a 
selection would include a wider variety of galaxies that 
are not so suitab le analogs of high r edshift objects com- 
pared to LBAs (jHoopes et al.|[2QQ7[ ). An example of a 
particularly clumpy LBA is the lum inous blue co mpact 
dwarf galaxy Haro 11 at z = 0.02 (Osthn et alJ|2009D, 
which lies close to the edge of our LBA selection cri- 
teria based on FUV luminosity and surface brightness 
(jGrimes et al.|[2QQ7l ). It is similar to the LBAs in most 
aspects but with a relatively high degree of dumpiness 
due to three strong light concentrations in an otherwise 
amorphous, forming galaxy. Its UV morphology at high 
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redshift is similar to that of double nucleated galaxies 
identified on t he basis of a lar^e M 20 or even a "by eye" 
classification (lOverzier et al.|[2QQ8[ ). 

4.3. Possible Implications for accretion processes in 
starbursts at high redshift 

Based on the results presented in the previous sections, 
we can now make a series of statements purely based on 
the comparison of the morphologies of LBAs and high 
redshift starbursts (such as BzKs and LBGs): 

1. LBAs show clumpy star formation. Our redshift 
simulations suggest that real high redshift samples have 
a similar (LBGs; z ~ 3 — 4) or perhaps a slightly stronger 
{BzKs; z ^ 2) degree of clumpy star formation. 

2. LBAs show luminous UV clumps and faint op- 
tical tidal features that, in the majority of cases, are 
interpreted by us as being due to mergers based on 
the visual evidence. Our simulations suggest that this 
information can not be recovered from the rest-frame 
UV morphological parameters at high redshift, while 
in the rest-frame optical perhaps 1 out of every 2 
objects would be classified as disturbed. This clearly 
demonstrates that current estimates for the fraction of 
galaxies at high redshift having disturbed morphologies 
(as opposed to smooth/symmetric morphologies) is most 
likely a (weak) lower limit on the true fraction. This 
statement is independent of the physical mechanism 
that is the cause of these disturbed morphologies. 

3. The observations suggest a scenario where either 
cold gas accretion and internal instabilities can drive 
clumpy star formation in LBGs at levels higher than 
any seen in LBAs in the local universe, or where clumpy 
cold accretion, possibly at the level of major or minor 
mergers may be responsible for star formation in LBAs 
and LBGs alike. We will discuss these possibilities in 
more detail below. 

The latest observational evidence suggests that the 
fraction of LBGs at z ~ 3 — 4 having dist urbed or dis- 
torted morp holo Ries is ^30% (e.g., Conseli ce et al.ll2QQ3l : 
iLotz et al.l [20061 : iConselice fc Arnold „2Q09. This Paper). 
Although these disturbances do not necessarily need to 
be explained by galaxy mergers, it does suggest that bulk 
material is coming in, p erhaps in the form o f giant gas 
clouds or minor mergers. iGonselice &: Arnoldl ([2009) find 
that the fraction of 2; ^ 4 5-dropouts forming a pair 
with another 5-dropout is very similar to the fraction 
of disturbed morphologies: ~20% within a 20 h~^ kpc 
projected radius. Interpreting the similarity in the pair 
counts and morphologies as evidence for merging, they 
estimate that a ~ 10^^ Mq galaxy will undergo a ma- 
jor merger every 1-2 Gyr at 2: ~ 3 — 4. Similar high 
merger rates were obtained for z ^ 2 — 3, possibly with 
an incre ase towards the most massive and most luniinous 
galaxies (i Consehc e et al.ll2003l : lCon sehce fc Arn"oldll2009l : 
IB luck et a l. 2009). With such high inferred merger rates, 
the mass growth of LBGs due to these mergers could 
make a significant contribution to the mass growth at 
z ~ 3 compared to that due to star formation fueled by 
a more continuous gas accretion and sustained over a pe- 
riod of a Gyr at a rate of ~ 10 Mq yr~^. Perhaps these 



two channels of formation are not contradictory given 
that the high density environments of LBGs are likely to 
be simultaneously associated with both frequent merg- 
ers of small galaxies and rapid cold accretion (likely with 

some level of lumpiness)^ 

TConselice & Arnold (i2009l ) suggest that the remain- 
ing 70-80% of LBGs, that appear as relatively 
smooth/symmetric and are forming stars at a similarly 
high rate as the disturbed objects, could be the result 
from rapid gas collapse, provided that a sufficiently long 
amount of time (^0.5 Gyr) has passed since their last ma- 
jor merger otherwise this would have been apparent in 
their morphology given the high (inferred) merger rates. 
However, our new results based on LBAs present an im- 
portant caveat to this interpretation: we have shown 
that, at least for galaxies at high redshift similar to 
LBAs, in the majority of cases we are not able to de- 
tect significant disturbances or asymmetries to their mor- 
phologies (mostly because of redshift effects, and not be- 
cause they are not there). This allows for the possibility 
that the fraction of galaxies at high redshift that is under- 
going interactions could be much higher than currently 
inferred from the observations. Alternatively, less violent 
accretion processes coupled with large disk instabilities 
at high redshift are capable of triggering starbursts at 
levels only seen in low redshift, merger-induced samples 
such as the LBAs. Future deep, high resolution observa- 
tions of rest-frame optical morphologies and kinematics 
(see below), together with constrai nts on the LBG pair 
fract i ons and small-scal e clustering (jConselice fc Arnoldl 
I2OO9I : iCooke et al.l 120091 ) will perhaps allow us to distin- 
guish between these scenarios. 

4.4. Relation to Studies of Gas kinematics 

In recent years, the study of high redshift LBGs has ad- 
vanced considerably beyond studies that are based purely 
on morphology. Mainly through the use of integral field 
spectrographs in the near-infrared it has become feasi- 
ble to determine the basic kinematical properties of the 
emission line gas as well. Motivated by the high de- 
gree of similarity in the properties of our locally selected 
UV-luminous galaxies and those at high redshift, it is 
a useful exercise to compare the gas ki nematics of LBAs 
and L BGs. In a first study published bv lBasu-Zych et aTl 
(|2009f ). the bright Pa-a emission line was used to study 
the resolved gas dynamics in 3 of our LBAs. In two cases, 
a mild velocity gradient was found, but in all three cases 
the kinematics were dominated by the dispersion rather 
than structured rotation {v/cf ^ 1 — 2). Simulating the 
data at z ~ 2 demonstrated that the (gas) kinematical 
properties of these objects are similar to the kinematical 
profiles commonly seen at high redshift (e.g. Law e t al.l 
2007, 2009; Genzel et al. 200^ iForster Schreiber eFall 
,,2009: Lowenthal et al. 200a). 

iLehnert et al.l (|2009l ) argued that neither the self- 
gravity of disks fueled by gas accretion fiows nor the 
internal velocity dispersions of massive star-forming 
clumps can fully explain the high velocity dispersions 
observed at high redshift. Furthermore, the most mas- 
sive of clumps observed possibly would not have been 
formed at all if the gas tu rbulence in the disk w as not 
high enough to begin with (jElmegreen et al.ll2009f ). while 
mergers alone may not be sufficient to explain those ob- 
jects dominated by a large number of massive clumps 
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in so-called ''clump-clust er" or ''chain galaxy" c onfi^u- 
rations at high redshift (|BQurnaud fc Elmegreen| [2QQ9). 
Instead, 'Lehnert et al." (2009') suggest that the kinematic 
properties of the ISM in starburst galaxies is affected by 
the mechanical energy input resulting from massive star 
formation. In this scenario, the starbursts that are as- 
sociated with each of the clumps drive blast waves from 
super novae (SN) and stellar winds that appear sufficient 
to give rise to the high gas pressures a nd hig;h veloc- 
ity dispersions observed. As shown in Ove rzier et al.l 
(f2009), some LBAs show very high pressures and clear 
evidence for an ISM that is dominated by starburst- 
and SN-feedback associated with massive star-forming 
clumps, consistent with such a scenario. If this is cor- 
rect, then estimates for the merger rates at high red- 
shift estimated from pair counts or morphologies are 
perhaps less biased than those obtained from the (ion- 
ized) gas kinematics because the latter may not always 
trac e galaxy interactions, if pres ent, very well. In addi- 
tion, [Robertson fc Bullockl (|2008 l) have shown with simu- 
lations that at least some merging systems would still be 
classified as "disks" using the methodology applied to ob- 
servatio ns at z ^ 2 by Shapiro et al. (2008). The simula- 
tions of lRobertson fc Bullockl (|2008 ) also show how some 
of the main structural, spectral and chemical properties 
of certain BzK galaxies can be explained in the con- 
text of gas-rich mergers as well, eventhough it has been 
claimed that such systems cannot be merger remnants. 
Similar studi es of the ga s kinematics in LBAs as initi- 
ated by Basu-Zvch et al.l ([2009) will be very useful, and 
will allow us to carefully test the methods typically ap- 
plied to high redshift starbursts in a suitable low redshift 
comparison sample. Results on the kinematical proper- 
ties in a much la rger sample of LBAs are forthcoming 
(jGoncalves et al.fe OlO. in prep.). 

5. SUMMARY 

• We have compared the sam ple of 30 nearj j y UV - 
selected starburst galaxies from i Qverzier et all ()2009i) 
with samples of starburst galaxies at z ^ 2, z ^ 3 and 
2: ~ 4 selected from the Hubble Ultra Deep Field ACS 
and WFC3 observations. 

• These so-called Lyman Break Analogs have compa- 
rable UV colors (a probe of dust), UV-optical colors (a 
probe of age), metallicity and half-light radii compared 
to BzK galaxies at 2; ~ 2 and LBGs at z ~ 3 — 4 all 
selected to have the same rest-frame UV luminosity (a 
probe of SFR) of Luv > 0.31'=^. 

• LBAs lie on a stellar mass-metallicity relation that 
is offset from typical local galaxies of the same mass, 
but similar to that observed for starburst galaxies at 
z ^ 2. This indicates that these starburst galaxies 
are still in the process of converting relatively pris- 
tine gas into stars or had a recent accretion event of 
metal-poor gas, possibly coupled with outflows of metals. 

• We have determined the morphological parameters 
(G, C, A, and M20) of LBAs in the rest-frame UV and 
optical, and performed redshift simulations to study 
the effects of degradations in physical resolution and 
sensitivity on morphological classifications. The high 
UV luminosities of the LBAs allow us to make such a 



comparison without the need for artificial brightening 
as employed in previous studies. While at low redshift 
most LBAs can be classified as being disturbed (only!) 
on the basis of a high asymmetry, there is a significant 
reduction in the asymmetries at high redshift. This 
reduction is less in the optical than in the UV. These 
results suggest that morphological disturbances in 
starburst galaxies similar to LBAs can be easily missed 
in current observations of high redshift galaxies. 

• We have compared the morphologies of BzKs 
and LBGs with those of LBAs simulated at a similar 
redshift. For the rest-frame UV and optical comparison, 
the LBAs were simulated in ACS/WFC Vqqq and 
WFC3/IR i^i6o images, respectively, having similar 
exposure times as the actual HUDF observations. The 
measured morphologies are generally very similar for 
the three samples with a few exceptions: the median 
M20 of LBGs at 2: ~ 3 — 4 is larger than that of LBAs, 
while BzK galaxies have the highest A and M20 and 
the lowest C and G, consistent with a higher degree of 
dumpiness. BzK galaxies are also somewhat redder and 
larger than LBGs (and LBAs), consistent with previous 
findings. 

• It has been suggested that high redshift galaxies 
experience intense bursts unlike anything seen in the 
local universe, possibly due to cold flows and instabili- 
ties. In part, this is based on the fact that the majority 
(~70%) of LBGs do not show morphological signatures 
of interactions or mergers. Our results suggest that this 
evidence is insufficient, since in the majority of cases 
we are not able to detect significant disturbances or 
asymmetries in LBAs artificially redshifted to z ~ 2 — 4. 
Likewise, some conclusions drawn from the nebular 
gas kinematics in LBGs have also been shown to be 
ambiguous. This leaves open the possibility that, at 
least in starburst galaxies such as the ones discussed 
here, clumpy accretion and mergers remain important 
processes, possibly together with rapid gas accretion 
through other means. 

As advocated in this paper, the HST sample of LBAs is 
extremely well-suited for performing morphological com- 
parisons with other samples of galaxies at a wide red- 
shift range, or, for example, for testing morphological 
classification schemes. In the next decade we can expect 
significant improvements in the measurements of merger 
rates, morphologies, gas fractions and gas and stellar 
kinematics in large samples of high redshift galaxies us- 
ing the James Webb Space Telescope (JWST), the At- 
acama Large Millimeter Array (ALMA), the Extremely 
Large Telescope (ELT), and the Square Kilometer Array 
(SKA). Also, improved semi-analytic modeling of high 
redshift galaxies and detailed hydro-simulations of the 
formation and evolution of clumpy systems will better 
constrain the importance of mergers, infall and general 
dissipational processes as a function of redshift. The 
study of similar processes occurring in nearby galaxies 
such as the LBAs studied in this paper provides an in- 
valuable tool for comparing with the observational record 
at high redshift. The authors would be happy to make 
the calibrated images available on request. 
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TABLE 2 

Summary of Medians and Scatters"^ for the Morphological Parameters from Figs. [SHU 









Rest-UV 






Rest- 


optical 




Sample 




G 


C 


M20 


A 


G 


C 


M20 


A 


LBA 


0.2 


0.60 ±0.04 


3.13 ±0.57 


-1.56 ±0.36 


0.44 ±0.15 


0.61 ±0.05 


3.62 ±1.38 


-1.56 ±0.77 


0.44 ±0.15 




2.0 


0.58 ±0.05 


3.15 ±0.51 


-1.74 ±0.23 


0.25 ±0.09 


0.54 ±0.04 


2.89 ±0.34 


-1.77 ±0.20 


0.36 ±0.10 




3.0 


0.59 ±0.04 


3.18 ±0.49 


-1.79 ±0.25 


0.22 ±0.07 


0.52 ±0.03 


2.90 ±0.40 


-1.76 ±0.26 


0.36 ±0.15 




4.0 


0.55 ±0.05 


3.09 ±0.46 


-1.83 ±0.29 


0.20 ±0.05 










BzK 


2.0 


0.54 ±0.11 


2.51 ±0.73 


-1.12 ±0.42 


0.34 ±0.13 


0.49 ± 0.08 


2.58 ±0.31 


-1.60 ±0.19 


0.29 ±0.18 


LBG 


3.0 


0.52 ±0.18 


2.69 ±0.99 


-1.36 ±0.55 


0.23 ±0.14 


0.50 ± 0.04 


2.72 ±0.27 


-1.65 ±0.22 


0.25 ±0.18 


LBG 


4.0 


0.52 ±0.18 


2.75 ± 1.05 


-1.44 ±0.59 


0.14 ±0.11 











' The standard deviation on the mean is estimated from the bias-corrected sample variance. Note, however, that the distributions are often not symmetric 
around the median or mean. 



TABLE 3 

Statistical Significance^ of the Differences between Morphological 
Parameters from Figs. EHZIin the rest-frame UV. 



Sample (z) 


G 

R KS 


C 

R KS 


M20 
R KS 


A 

R KS 


LBA-LBA 0.2^2.0 
0.2^3.0 
0.2^4.0 
LBA-BzK 2.0^2.0 
LBA-LBG 3.0^3.0 
LBA-LBG 4.0^4.0 


-0.34 3e-l 
-0.24 le-1 
-0.77 le-2 
-0.35 5e-3 
-0.39 le-3 
-0.17 8e-2 


0.02 5e-l 
0.07 2e-l 
-0.05 5e-l 
-0.72 6e-4 
-0.45 3e-4 
-0.29 le-2 


-0.43 5e-2 
-0.51 5e-2 
-0.58 le-1 
1.30 6e-8 
0.71 3e-6 
0.59 7e-4 


-1.10 5e-6 
-1.32 2e-9 
-1.54 6e-ll 
0.60 le-2 
0.08 2e-l 
-0.52 3e-4 



' We define the significance using the parameters R and KS. R = (P2 — ^l) / y '^^^ "l~ '^^^^2 ' ^^^''^ 
and cTp. are the median and standard deviation of the morphological parameter P measured for sample 
i. The standard deviation is estimated from the bias-corrected sample variance. The parameters KS 
indicate the value of significance that the null hypothesis is true for the two-sided K-S statistic applied 
to each pair of samples. 



TABLE 4 

Statistical SiGNiFicANCEt of the Differences between Morphological 
Parameters from Figs. [SHZIin the rest-frame optical. 



Sample (z) 


G 

R KS 


C 

R KS 


M20 
R KS 


A 

R KS 


LBA-LBA 0.2^2.0 
0.2^3.0 
0.2^4.0 
LBA-BzK 2.0^2.0 
LBA-LBG 3.0^3.0 
LBA-LBG 4.0^4.0 


-1.10 2e-5 
-1.49 6e-8 

-0.58 6e-5 
-0.35 5e-2 


-0.52 2e-3 
-0.50 2e-3 

-0.66 3e-4 
-0.38 4e-2 


-0.26 5e-3 
-0.25 le-1 

0.61 7e-4 
0.32 5e-2 


-0.46 5e-2 
-0.39 le-1 

-0.36 3e-l 
-0.47 5e-2 



sing the pa 



and CTp. are the median and standard deviation of the morphological parameter P measured for sample 
i. The standard deviation is estimated from the bias-corrected sample variance. The parameters KS 
indicate the value of significance that the null hypothesis is true for the two-sided K-S statistic applied 
to each pair of samples. 
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